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HIGHLIGHTS 


GRAPHICAL  ABSTRACT 


•  We  characterize  layered  MoS2  as  an 
anode  for  sodium-ion  batteries. 

•  The  strong  binding  of  Na  with  MoS2 
is  favorable  against  Na  cluster 
formation. 

•  Sodiation  induces  a  phase  trans¬ 
formation  from  hexagonal  to  tetrag¬ 
onal  phase. 

•  We  found  the  maximum  capacity  of 
146  mAh  g~]  and  average  electrode 
potential  of  0.75-1.25  V  for  Na 
intercalation  in  MoS2. 

•  The  result  on  Na  diffusion  barrier  in 
MoS2  suggests  a  moderately  fast 
mobility. 
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Identifying  suitable  layered  materials  as  electrodes  with  desirable  electrochemical  properties  remains  a 
key  challenge  for  rechargeable  Na-ion  batteries  (NIBs).  Using  first  principles  methods,  here  we  examine 
the  efficacy  of  layered  molybdenum  disulphide  (MoS2)  as  a  host  electrode  material  for  NIBs.  We  identify 
various  low  energy  Na  adsorption  sites  and  evaluate  the  stability  of  the  hexagonal  and  tetragonal  pol¬ 
ytypes  of  MoS2  upon  Na  intercalation.  Our  results  illustrate  a  moderately  strong  binding  between  Na  and 
MoS2  that  is  thermodynamically  favorable  against  the  cluster  formation  and  phase  separation  of  Na.  We 
find  that  while  Na  intercalation  in  MoS2  results  in  a  phase  transformation  from  the  hexagonal  phase  to 
the  tetragonal  phase,  it  gives  rise  to  a  maximum  theoretical  capacity  of  146  mAh  g_1  and  a  low  average 
electrode  potential  in  the  range  of  0.75-1.25  V.  Our  calculations  of  Na  diffusion  kinetics  indicates  a 
moderately  fast  mobility  of  Na  in  the  van  der  Waals  interlayer  spaces  of  MoS2.  These  results  highlight  the 
promise  of  MoS2  as  an  appealing  negative  electrode  (anode)  for  rechargeable  NIBs. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 
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Rechargeable  Na-ion  batteries  (NIBs)  are  rapidly  emerging  as  a 
plausible  alternative  to  Li-ion  batteries  (LIBs),  particularly  for  sta¬ 
tionary  grid-based  energy  storage  [1—3],  This  growing  interest  in 
NIBs  can  be  attributed  not  only  to  the  abundance  of  Na  and  the 
associated  cost  advantages,  but  also  to  its  low  electronegativity 
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which  is  next  only  to  Li  among  alkali  metals.  An  increasing  number 
of  experimental  and  theoretical  studies  have  investigated  a  wide 
range  of  Na-ion  battery  chemistries  in  pursuit  of  desirable  elec¬ 
trochemical  capacities,  an  appropriate  voltage  range  and  a  long 
cycle  life  [1,2].  A  significant  part  of  this  effort  is  focused  on  iden¬ 
tifying  suitable  van  der  Waals-bonded  layered  materials  that  can 
act  as  negative  electrodes  (anodes)  of  NIBs  [2],  In  contrast  to  the 
potential  materials  such  as  Sn  and  Pb  which  offer  large  capacities 
but  poor  mechanical  stability  [4],  layered  materials  can  intrinsically 
accommodate  guest  ions  within  the  van  der  Waals  spacing  with  a 
minimum  structural  distortion  and  allow  for  their  fast  transport 
during  charging  and  discharging  cycles.  In  state  of  the  art  LIBs, 
graphitic  carbon  is  the  most  commonly  used  layered  material  for 
anodes  [2  .  Li  forms  an  intercalation  compound  LiC6  with  graphite 
at  full  charge  with  a  maximum  capacity  in  excess  of  340  mAh  g_1 
[5],  However,  due  to  the  larger  size  of  Na  ion  (1.02  A  vs.  0.76  A  for 
Li),  Na  intercalation  in  graphite  is  highly  inefficient — Na  can  only 
form  an  eighth-stage  intercalation  compound  NaC64  with  graphite 
with  a  very  poor  capacity  [6], 

Among  the  non-carbon  layered  materials  that  are  being 
explored  as  effective  electrode  materials  for  NIBs,  titanates  and  ti¬ 
tanium  phosphates  have  demonstrated  a  maximum  capacity  in  the 
range  of  100—140  mAh  g_1  [2,3  .  However,  a  large  class  of  layered 
materials,  namely,  transition  metal  dichalcogenides  remains  rela¬ 
tively  unexplored  as  a  potentially  useful  alternative  to  graphite  for 
anodes  of  NIBs.  These  materials  (with  chemical  formula  MX2, 
where  M  is  a  group  IV,  V  or  VI  transition  metal  such  as  Mo,  W  and  Ti, 
and  X  is  a  chalcogen  such  as  S,  Se  and  Te)  are  of  interest  for  a  broad 
range  of  applications  in  hydrogen  storage,  catalysis,  and  more 
recently,  in  optoelectronic  devices  [7,8],  Similar  to  graphite,  the 
atomic  structure  of  this  family  of  materials  is  characterized  with  a 
strong  in-plane  and  a  weak  out-of-plane  bonding  9,10  .  Moreover, 
transition  metal  dichalcogenides  have  a  larger  interlayer  spacing 
than  graphite,  and  offer  a  wide  range  of  electron  affinity  for  the 
intercalation  of  a  guest  electronegative  ion  [11], 

While  this  ability  to  readily  intercalate  guest  ions  has  been 
greatly  exploited  for  exfoliating  these  layered  structures  into  indi¬ 
vidual  layers  [12],  attempts  to  employ  transition  metal  dichalco¬ 
genides  as  active  electrode  materials  for  rechargeable  batteries 
have  been  few,  and  so  far,  largely  limited  to  LIBs  [13—15].  Recently, 
Park  et  al.  demonstrated  electrochemical  insertion  of  Na  in  layered 
molybdenum  disulphide  (M0S2)  with  a  steady  state  capacity  of 
85  mAh  g_1  after  100  charge— discharge  cycles  [16],  Similarly, 
studies  have  also  reported  a  specific  capacity  in  the  range  of 
50—140  mAh  g_1  for  Nb  [17],  W  [18],  and  Ti  [19]  disulphides.  While 
these  initial  studies  report  moderate  capacities  of  transition  metal 
dichalcogenides  for  Na,  very  recent  experiments  have  shown  that 
they  can  be  further  enhanced  to  250—350  mAh  g_1  by  using 
composites.  For  instance,  the  work  of  David  et  al.  [20]  and  Bang 
et  al.  [21]  demonstrated  that  when  a  composite  of  M0S2  with 
graphene  oxide  is  employed  as  an  active  anode,  the  maximum 
capacity  of  230  and  360  mAh  g  ,  respectively,  could  be  achieved. 
These  experiments  point  to  the  feasibility  of  layered  transition 
metal  dichalcogenides  as  high  capacity  anode  materials  for  NIBs, 
however,  an  atomic  level  understanding  of  Na  intercalation  pro¬ 
cesses  in  these  materials  is  not  yet  available.  The  knowledge  of 
structural  stability  of  Na-intercalated  transition  metal  dichalcoge¬ 
nides,  their  electronic  properties  and  electrochemical  characteris¬ 
tics  over  a  wide  compositional  range  is  essential  in  determining 
their  viability  as  electrode  materials  for  NIBs. 

In  this  work,  we  systematically  investigate  Na  intercalation  in 
layered  M0S2  as  a  representative  example  of  layered  transition 
metal  dichalcogenide  anode  for  NIBs.  We  consider  all  possible 
structural  polytypes  that  have  been  typically  reported  for  crystal¬ 
line  M0S2.  In  its  natural  state,  M0S2  is  generally  found  in  its  most 


stable  form  with  a  hexagonal  2H  crystal  structure  [11],  while  a  high 
energy  tetragonal  IT  polytype  can  also  be  stabilized  via  liquid-state 
exfoliation  [22].  Moreover,  several  experimental  and  theoretical 
studies  have  pointed  out  a  coexistence  of  2H  crystal  symmetry  with 
a  closely  related  rhombohedral  3R  symmetry  [23,24],  while  a  few 
studies  have  also  investigated  a  hypothetical  tetragonal  2T  struc¬ 
ture  [25],  Using  first  principles  density  functional  theory  calcula¬ 
tions,  here  we  analyze  the  structural  stability,  electronic  structure, 
and  intercalation  voltages  for  Na  intercalation  compounds  with  2H, 
3R,  as  well  as  IT  and  2T  polytypes  of  M0S2.  We  find  that  Na  inter¬ 
calation  in  M0S2  leads  to  a  structural  phase  transformation  from 
the  hexagonal  to  the  tetragonal  polytype,  and  further,  all  interca¬ 
lated  compounds  remain  largely  metallic  irrespective  of  the 
structural  polytype.  We  obtain  crucial  electrochemical  properties 
such  as  average  operating  voltage,  theoretical  capacity  and  Na 
diffusion  barriers,  and  evaluate  the  feasibility  of  layered  M0S2  as  an 
effective  anode  material  for  NIBs. 

2.  Methodology 

Schematic  representations  of  the  crystal  structures  for  the  2H, 
3R,  IT  and  2T  MoS2  polytypes  considered  in  our  study  are  shown  in 
Fig.  1.  Both  2H  and  3R  polytypes  are  characterized  by  a  trigonal 
prismatic  D3/,  coordination,  while  IT  and  2T  polytypes  have  an 
octahedral  Of,  coordination  [11],  MoS2  layers  in  2H,  3R,  IT  and  2T 
polytypes  are  arranged  in  antiparallel  AB,  parallel  ABC,  parallel  AA, 
and  antiparallel  AA'  stacking  sequence,  respectively.  For  all  these 
ciystal  symmetries,  we  performed  first  principles  total  energy 
calculations  within  the  framework  of  spin-polarized  density  func¬ 
tional  theory  (DFT)  using  Vienna  Ab  initio  Simulation  Package 
(VASP)  [26].  We  used  the  generalized  gradient  approximation  of 
the  Perdew— Burke— Ernzerhof  form  to  describe  electron  exchange 
and  correlation  [27],  while  the  core  and  valence  electronic  in¬ 
teractions  were  modeled  using  projector  augmented  wave  poten¬ 
tials  [28],  The  electrons  in  2p63s1,  4p65s14d5,  2s2 3p4  states  were 
treated  as  valence  electrons  for  Na,  Mo  and  S,  respectively.  Since  the 
standard  generalized  gradient  approximation  underestimates 
interlayer  van  der  Waals  interactions,  we  employed  Grimme's  en¬ 
ergy  correction  method  to  correctly  account  for  these  dispersion 
interactions  [29  .  A  17  x  17  x  11  fc-point  mesh  per  M0S2  unit  for 
sampling  the  Brillouin  zones  was  used  to  ensure  the  variation  in 
total  energy  is  less  than  1  meV/atom.  In  each  case,  we  used  a  plane 
wave  kinetic  energy  cut  off  of  500  eV.  The  optimized  structures 
were  obtained  by  allowing  both  the  cell  vectors  as  well  as  ionic 
positions  to  fully  relax  until  the  Hellmann-Feynman  forces  were 
less  than  0.01  eV  A-1. 

Our  calculated  optimized  lattice  parameters  for  all  M0S2  poly¬ 
types  are  tabulated  in  Table  SI,  SD.  The  obtained  lattice  parameters 
for  all  polytypes  are  in  good  agreement  with  experimental  and 
theoretical  values  reported  earlier  [11,25,30],  Consistent  with 
experimental  observations,  we  found  that  2H  is  energetically  the 
most  favorable  polytype  of  M0S2  (see  Table  SI  in  SD)  [24],  Finally, 
we  also  computed  total  density  of  states  for  all  polytypes  of  M0S2 
and  found  that  2H  and  3R  polytypes  are  semiconductors  with 
0.9  eV  and  0.92  eV  band  gap,  respectively,  while  IT  and  2T  poly¬ 
types  are  metallic  (Fig.  SI  in  SD).  These  observations  are  also 
consistent  with  earlier  reports  investigating  the  electronic  struc¬ 
ture  of  transition  metal  dichalcogenides  [31,32],  To  obtain  opti¬ 
mized  geometries  and  energies  of  Na  intercalation  compounds 
with  M0S2,  Na  atoms  were  first  inserted  in  the  interlayer  spacings 
to  form  the  compound  Na„MoS2,  and  the  composite  Na- 
intercalated  structure  was  reoptimized.  The  binding  energy  for  Na 
in  layered  M0S2  was  then  calculated  using  the  expression 
£b(x)  =  £(NaxMoS2)  -  £(MoS2)  -  x£(Na),  where  £(NaxMoS2)  is  the 
total  energy  of  the  Na*MoS2  compound,  £(MoS2)  is  the  total  energy 
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Fig.  1.  Structural  polytypes  of  layered  MoS2:  (a)  2H,  (b)  3R,  (c)  IT  and  (d)  2T.  a,  b  and  c  denote  the  lattice  constants  of  the  hexagonal  unit  cell.  Yellow  and  purple  spheres  represent  S 
and  Mo  atoms,  respectively.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 


of  the  same  M0S2  polytype  and  £(Na)  is  the  energy  of  Na  atom  in  its 
gaseous  state.  According  to  this  definition,  a  negative  binding  en¬ 
ergy  indicates  an  exothermic  chemical  interaction  between  Na  and 
M0S2. 

3.  Results  and  discussion 

3.3.  Na  adsorption  sites  and  energetics  of  intercalation  in  M0S2 

To  systematically  study  Na  intercalation  in  layered  M0S2,  we 
first  examined  low  energy  adsorption  or  binding  sites  for  an  iso¬ 
lated  Na  atom  in  the  interlayer  van  der  Waals  spacing.  We  identi¬ 
fied  two  distinct  Na  binding  sites  irrespective  of  the  M0S2  polytype, 
namely,  (a)  an  octahedral  site,  where  Na  is  bonded  to  six  nearest 
sulfur  atoms,  (b)  a  tetrahedral  site,  where  adsorbed  Na  is  bonded  to 
four  nearest  sulfur  atoms.  These  binding  sites  are  schematically 
illustrated  in  Fig.  2.  We  find  that  among  these  sites,  the  adsorption 
of  an  isolated  Na  is  more  favorable  in  octahedral  site— Na  binding 
energy  at  the  octahedral  site  is  generally  lower  by  80—300  meV 
than  at  the  tetrahedral  site. 

Next  we  investigated  the  concentration  dependence  of  Na 
intercalation  in  all  four  polytypes  of  M0S2  using  a  fully  periodic 


(a)  (b)  (c) 


Fig.  2.  Schematic  illustration  of  three  distinct  binding  sites  for  Na  in  the  interlayer 
spacing  of  MoS2:  (a)  octahedral,  (b)  tetrahedral,  and  (c)  distorted-octahedral.  Orange 
spheres  represent  Na  atoms.  (For  interpretation  of  the  references  to  color  in  this  figure 
legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 


supercell  consisting  of  at  least  two  M0S2  layers,  each  containing 
four  M0S2  units.  Using  such  a  large  supercell  allows  us  to  consider 
NaxMoS2  compounds  with  eight  distinct  chemical  compositions 
(x  =  0.125,  0.25,  0.375,  0.5,  0.625,  0.75,  0.875  and  1).  Since  Na 
ordering  can  be  different  for  different  Na  concentrations,  we 
considered  all  possible  configurations  for  every  composition 
considered.  In  all  configurations,  Na  atoms  were  initially  distrib¬ 
uted  in  either  octahedral  or  tetrahedral  sites  as  far  away  from  each 
other  as  possible  in  order  to  minimize  their  self-interaction.  For 
each  configuration,  the  structure  was  then  fully  optimized.  The  Na 
ordering  that  gives  the  lowest  binding  energy  for  each  concentra¬ 
tion,  the  ground  state,  was  chosen  for  further  electrochemical 
analysis.  Our  calculations  show  that  both  geometrical  characteris¬ 
tics  and  the  energetics  of  Na-intercalated  3R  and  2T  M0S2  are 
almost  identical  to  the  geometry  and  the  energetics  of  2H  and  IT 
polytype,  respectively.  This  observation  suggests  that  the  weak 
dispersion  interactions  that  distinguish  2H  and  3R  polytypes,  as 
well  as  IT  and  2T  polytypes,  are  negligible  compared  to  the  strong 
chemical  interactions  between  Na  and  M0S2.  Therefore,  here  we 
present  the  results  for  Na  intercalation  in  2H-  and  IT-M0S2  as 
representative  M0S2  polytypes,  while  key  the  results  for  3R  and  2T 
M0S2  polytypes  are  summarized  in  the  SD. 

Fig.  3  shows  Na  binding  energy  across  the  entire  composition 
range  for  various  stable  binding  configurations  within  2H-  and  1T- 
M0S2.  It  also  presents  a  convex  hull  for  the  binding  energies  (black 
lines),  representing  lowest  energy  binding  sites  in  both  2H  and  IT 
polytypes.  Interestingly,  in  addition  to  the  octahedral  and  tetrahe¬ 
dral  binding  sites,  we  also  identified  a  distinct  distorted-octahedral 
Na  binding  site  in  several  configurations.  The  distorted-octahedral 
binding  site  is  structurally  close  to  the  octahedral  site  but  with  a 
Na— S  polyhedron  that  is  distorted  from  its  original  octahedral 
shape  due  to  a  relative  sliding  of  S  atomic  plane  (see  Fig.  2c).  It  is 
evident  that  Na  binding  energies  for  all  configurations  in  2H-  and 
IT-M0S2  polytypes  are  negative,  indicating  an  energetically  favor¬ 
able  chemical  interaction  between  Na  and  layered  M0S2.  At  low 
(x  <  0.25)  and  full  (x  =  1)  Na  concentrations  in  2H-NaxMoS2, 
octahedral  sites  provide  lowest  energy  configurations,  while 
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tetrahedral  sites  correspond  to  the  metastable  configurations.  In 
contrast,  for  intermediate  concentrations  in  2H-MoS2,  distorted- 
octahedral  sites  provide  the  energetically  most  favorable  Na  bind¬ 
ing  configurations.  On  the  other  hand  for  nearly  all  Na  concentra¬ 
tions  in  1T-MoS2,  distorted-octahedral  sites  are  found  to  be 
energetically  most  favorable.  Fig.  3  in  SD  shows  the  lowest  energy 
configurations  of  both  2H-  and  lT-Na*MoS2,  while  their  optimized 
lattice  parameters  are  summarized  in  Table  1. 

While  Na  binding  energies  in  2H-  and  IT-M0S2  are  negative  for 
all  concentrations  as  shown  in  Fig.  3,  it  is  clear  that  Na  binds  more 
strongly  to  IT-M0S2  than  to  2H-M0S2 — Na  binding  energies  in  1T- 
M0S2  are  generally  lower  by  1.1— 1.9  eV  than  the  binding  energies  in 
2H-M0S2.  The  stronger  Na  binding  with  the  IT  polytype  can  be 
explained  by  the  crystal  field  theory  that  describes  a  symmetry- 
induced  splitting  of  Mo-4d  orbitals  in  various  M0S2  polytypes 
[33,34],  In  2H-M0S2,  since  the  lowest  energy  Mo  -  4 df  orbital  of 
the  D3/l-MoS6  unit  is  fully  occupied,  an  electron  from  donor  Na  can 
only  be  accommodated  in  the  higher  energy  4dxy  and  4dx2_y2  or¬ 
bitals.  In  contrast,  the  lowest  energy  degenerate  4 dxy,  4 dyz  and  4dxz 
orbitals  in  IT-M0S2  are  partially  occupied  and  can  readily  accom¬ 
modate  the  electron  from  Na,  leading  to  a  large  negative  binding 
energy.  Nevertheless,  it  should  be  emphasized  that  Na  binding 
energies  for  all  concentrations  in  both  2H  and  IT-M0S2  polytypes 
are  much  lower  than  the  cohesive  energy  of  Na  (-1.11  eV  per  atom). 
This  observation  suggests  a  strong  driving  force  for  the  adsorption 
of  Na  within  the  interlayer  spaces  of  M0S2  and  against  the  forma¬ 
tion  of  Na  clusters  which  can  be  detrimental  to  the  performance  of 
electrochemical  cells.  Our  obtained  Na  binding  energies  therefore 
point  to  a  thermodynamic  suitability  of  layered  M0S2  as  an  elec¬ 
trode  material  for  NIBs. 

3.2.  Lattice  deformation  upon  Na  intercalation 

It  has  been  established  that  repetitive  structural  changes  asso¬ 
ciated  with  ion  intercalation  and  deintercalation  can  lead  to  ca¬ 
pacity  fade  and  a  loss  of  electrochemical  performance  of  the 
electrodes  [35],  To  investigate  the  influence  Na  intercalation  on  the 
stability  of  the  crystalline  structure  of  layered  M0S2,  we  have 
analyzed  the  variation  of  lattice  parameters  with  Na  concentration 
as  shown  in  Fig.  4.  We  find  that  for  both  hexagonal  and  tetragonal 
polytypes,  the  lattice  parameter  a  increases  only  marginally  in  the 
range  of  6.4— 6.8  A  over  the  entire  concentration  range.  In  contrast, 
the  effect  of  Na  intercalation  on  the  lattice  parameter  c  is  more 
significant.  As  Na  enters  the  van  der  Waals  interlayer  spaces,  c 
rapidly  rises  during  the  second  (x  =  0.125)  and  the  first  (x  =  0.25) 
stage  of  Na  intercalation,  and  fluctuates  in  a  narrow  range  upon 


Table  1 

Optimized  lattice  parameters  (a.  b,  c,  a,  (1,  7),  minimum  in-plane  Na-Na  distance 
(dNa— Na),  formation  energy  (Ef)  and  net  charge  on  Na  (q)  for  Na-intercalated  MoS2. 


X 

a.b  (A) 

c(A) 

a  (°) 

/?(“) 

y  (°) 

^Na-Na  (A) 

E/(eV) 

?(e) 

2H-NaxMoS2 

0.000  6.39 

12.40 

90 

90 

60 

0.125 

6.42 

13.58 

90 

90 

60 

6.42 

-0.11 

0.86 

0.250 

6.44 

14.87 

90 

90 

60 

6.44 

-0.21 

0.85 

0.375 

6.43 

15.27 

78 

77 

60 

5.46 

-0.33 

0.84 

0.500 

6.44 

14.84 

90 

90 

60 

3.19 

-0.47 

0.85 

0.625 

6.51 

15.10 

90 

102 

60 

3.25 

-0.56 

0.83 

0.750 

6.54 

14.68 

90 

90 

60 

3.20 

-0.66 

0.83 

0.875 

6.60 

14.98 

102 

89 

60 

3.21 

-0.72 

0.82 

1.000 

6.58 

14.87 

90 

90 

60 

3.29 

-0.73 

0.78 

lT-NaxMoS2 

0.000  6.46 

12.06 

90 

90 

60 

0.125 

6.52 

13.56 

90 

90 

60 

6.52 

0.38 

0.87 

0.250 

6.57 

14.35 

90 

96 

60 

6.50 

0.01 

0.87 

0.375 

6.58 

14.10 

88 

91 

60 

3.18 

-0.26 

0.86 

0.500 

6.60 

14.26 

89 

95 

60 

3.23 

-0.53 

0.85 

0.625 

6.67 

14.28 

88 

88 

59 

3.31 

-0.73 

0.84 

0.750 

6.74 

14.49 

76 

76 

60 

2.99 

-0.92 

0.84 

0.875 

6.77 

14.06 

89 

89 

59 

3.21 

-1.17 

0.83 

1.000 

6.82 

14.33 

103 

89 

60 

3.37 

-0.83 

0.82 

further  increasing  Na  concentration.  These  small  fluctuations  are 
associated  with  both  interlayer  and  intralayer  the  Na— Na  in¬ 
teractions.  As  a  consequence  of  interlayer  expansion,  Na  interca¬ 
lation  leads  to  a  volume  expansion  of  27%  at  full  concentration.  It 
should  be  noted  that  while  this  magnitude  of  volume  expansion  is 
apparently  large,  Na  intercalation  in  M0S2  in  the  concentration 
range  0  <  x  <  1  is  mainly  contributed  by  the  expansion  of  interlayer 
van  der  Waals  space,  and  not  due  to  any  bond  breaking  in  the  host 
material.  This  behavior  is  in  contrast  with  intercalation  in  the  bulk 
high  capacity  materials  for  rechargeable  batteries  (such  as  silicon 
and  tin),  where  the  intercalation  of  the  guest  ion  induces  bond 
breaking  and  a  subsequent  capacity  loss  in  the  host  material 
[36-39], 

3.3.  Phase  transformation  in  layered  MoS2  upon  Na  intercalation 

Next  we  evaluate  the  relative  stability  of  various  M0S2  polytypes 
upon  Na  intercalation.  It  is  well  known  that  among  all  M0S2  poly¬ 
types,  2H-M0S2  is  the  most  stable  form  of  M0S2.  We  therefore  define 
the  formation  energy  for  all  NaxMoS2  intercalation  compounds  with 
respect  to  2H-M0S2  and  Na  body  centered  cubic  crystal  as  the 
reference  states,  as  E/(x)  =  £(NaxMoS2)  -  £(2H  -  M0S2)  -  x£(Nabcc), 
where  £(NaxMoS2)  is  the  energy  of  NaxMoS2  compound  per  M0S2 


Fig.  3.  The  convex  hull  of  Na  binding  energies  for  (a)  2H-  and  (b)  lT-NaxMoS2.  Distinct  binding  energies  for  the  same  chemical  composition  of  NaxMoS2  indicate  metastable  Na 
binding  configurations. 
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x  in  NaxMoS2 


Fig.  4.  Variation  of  the  lattice  parameters  and  the  cell  volume  with  Na  concentration  for  2H-  and  lT-Na,MoS2  intercalation  compounds. 


formula  unit,  £(2H-MoS2)  is  the  energy  of  2H-M0S2  per  M0S2  for¬ 
mula  unit,  and  E(Nabcc)  is  the  energy  per  atom  of  Na  in  the  body 
centered  cubic  crystal.  A  negative  value  of  the  formation  energy 
indicates  a  thermodynamic  stability  of  the  intercalation  compound 
against  the  separation  into  most  stable  forms  of  both  M0S2  and  Na. 
The  calculated  values  of  the  formation  energy  for  the  lowest  energy 
configurations  for  2H-NaxMoS2  and  lT-NaxMoS2  are  presented  in 
Table  1,  while  Fig.  5  shows  the  variation  of  the  formation  energy  per 
Na  as  a  function  of  the  concentration  x.  It  can  be  seen  that  2H-M0S2  is 
stable  upon  Na  intercalation  over  the  entire  concentration  range, 
while  IT  polytype  is  stable  only  beyond  concentration  x  >  0.25.  More 
interestingly,  for  Na  concentrations  larger  than  x  =  0.39,  the  for¬ 
mation  energy  of  lT-NaxMoS2  is  lower  than  that  of  2H-NaxMoS2, 
indicating  a  phase  transformation  from  2H  to  IT  polytype  atx  =  0.39. 
This  phase  transformation  is  qualitatively  similar  to  the  observations 
reported  in  earlier  studies  on  the  intercalation-induced  phase 
transformation  in  transition  metal  dichalcogenide  nanotubes  and 
multilayers  [33,34,40].  These  studies  showed  that  the  phase  trans¬ 
formation  can  be  attributed  to  a  modification  of  Mo-d  orbitals  that 
stabilizes  octahedral  Mo  coordination  of  the  IT  phase  at  large  con¬ 
centrations  of  the  intercalating  electron  donors.  Recent  real-time 
high  resolution  microscopy  experiments  have  revealed  the  mecha¬ 
nism  that  such  phase  transformations,  which  involves  gliding  of  the 
S  atomic  plane  in  two  stages,  namely,  nucleation  of  IT  phase  and 
growth  of  IT  phase  via  migration  of  two  different  domain  bound¬ 
aries  [41,42]. 


Fig.  5.  Formation  energy  per  Na  atom  for  2H-  and  lT-NaxMoS2  as  a  function  of  Na 
concentration. 


3.4.  Charge  transfer  and  electronic  structure  of  NaxMoS2 

Next  we  analyze  the  nature  of  chemical  bonding  between  Na 
and  host  M0S2  and  the  electronic  structure  of  the  intercalation 
compounds.  Fig.  S4  in  the  SD  shows  a  typical  distribution  of  the 
charge  density  difference  in  Na  intercalated  layered  M0S2.  From 
the  charge  density  distribution,  the  net  charge  on  Na  was  obtained 
using  Bader  charge  population  analysis.  Table  1  shows  the  net 
charge  on  Na  for  low-energy  configurations  of  2H-  and  1T-Nax_ 
MoS2.  In  each  case,  it  is  evident  that  Na  donates  a  large  fraction  of 
its  valence  electron  indicating  a  strong  ionic  character  of  the  bond 
between  Na  and  the  host  M0S2.  Next,  we  calculated  the  total 
density  of  states  to  analyze  the  changes  in  electronic  structure  of 
M0S2  upon  Na  intercalation  as  a  result  of  the  charge  transfer.  Fig.  6 
shows  the  concentration-dependent  density  of  states  for  2H-  and 
lT-NaxMoS2.  It  is  clear  that  the  electronic  structure  of  2H-NaxMoS2 
switches  from  semiconducting  to  metallic  upon  initial  Na  inter¬ 
calation,  and  thereafter  remains  metallic  upon  further  intercala¬ 
tion.  On  the  other  hand,  the  metallic  IT-M0S2  is  found  to  remain 
metallic  over  a  large  concentration  range.  Interestingly,  at  full 
concentration  (x  =  1)  IT-M0S2  demonstrates  a  band  gap  in  excess 
of  0.98  eV.  This  unexpected  phenomena  can  be  attributed  to  the 
distortion  of  Mo  coordination  at  very  large  concentration  of 
available  electrons.  Similar  unexpected  changes  in  electronic 
structures  during  intercalation  have  also  been  observed  in  a  few 
electrode  materials  for  LIBs,  and  can  potentially  be  mitigated  by 
using  dopants  to  increase  conductivity  of  the  host  material 
[43-45], 

3.5.  Electrochemical  properties  of  Na-intercalated  M0S2 

In  order  to  examine  the  suitability  of  layered  M0S2  as  an  elec¬ 
trode  material  for  NIBs,  we  next  evaluate  its  crucial  electrochemical 
characteristics  such  as  voltage  profile,  specific  capacity  and  Na 
diffusion  barriers.  For  every  concentration  x  of  the  NaxMoS2  com¬ 
pound,  the  electrode  potential  V  with  respect  to  Na/Na+  is  calcu¬ 
lated  as  [46,47], 

E(x2)-E(x,)-(x2-x1)E(Nabcc) 

V  = - 5 - L 

e(x2  -xj) 

where  E(x2)  and  E(xi)  are  the  total  energies  of  the  NaxMoS2  com¬ 
pound  at  two  adjacent  low-energy  concentrations  x2  and  Xi  along 
the  computed  convex  hull  shown  in  Fig.  3,  and  E(Nabcc)  is  the  en¬ 
ergy  per  atom  of  Na  in  the  bcc  crystal.  Fig.  7  shows  the 
concentration-dependent  profile  of  the  electrode  potential  for  both 
2H  and  IT-M0S2  polytypes.  We  find  that  the  electrode  potential  for 
2H-NaxMoS2  varies  in  a  range  of  0—1  V  with  a  decreasing  trend  with 
increasing  capacity,  with  an  average  value  of  0.75  V.  This  observa¬ 
tion  is  qualitatively  similar  to  the  measurements  of  electrode 
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Energy  (eV) 


Fig.  6.  Total  density  of  states  for  (a)  2H-Na.MoS2  and  (b)  lT-Na2MoS2.  The  Fermi  level  is  set  to  zero  eV  in  each  case. 


potential  upon  discharge  as  reported  in  the  recent  work  on  Na 
intercalation  in  M0S2  [16],  On  the  other  hand,  the  voltage  profile  for 
lT-NaxIVIoS2  as  shown  in  Fig.  7(b)  varies  in  a  broader  range.  The 
larger  magnitude  of  the  concentration-dependent  voltages  for  1T- 
NaxMoS2  can  be  ascribed  to  the  stronger  binding  of  Na  with  1T- 
M0S2.  However,  as  shown  in  Fig.  5,  lT-NaxMoS2  compounds  are 
unstable  for  concentrations  x  <  0.25.  Therefore,  the  large  magni¬ 
tudes  of  the  potential  of  lT-NaxMoS2  for  low  concentrations  are 
likely  to  remain  inaccessible.  Considering  stable  compositions  of 
lT-NaxMoS2,  the  average  potential  is  obtained  as  1.25  V.  Interest¬ 
ingly,  the  average  potential  of  Na-intercalated  M0S2  is  much  lower 
than  Li-intercalated  M0S2  (2.0  V)  due  to  the  relatively  weaker 
binding  of  Na  compared  with  Li  [15,25],  Since  ideally  a  good  anode 
should  have  a  low  electrode  potential,  our  calculated  voltage  profile 
in  Fig.  7  suggests  that  layered  M0S2  is  suitable  as  an  anode  for  an 
NIB.  When  this  Na-intercalated  M0S2  anode  is  combined  with  high 
capacity  cathode  materials  such  as  Nai.5VPC>4.8Fo.7  and  Na3Mn- 
PO4CO3,  the  Na-ion  battery  cell  can  yield  a  desirable  open  circuit 
voltage  in  the  range  of  2.5— 3.5  V  [48,49]. 

From  our  calculations,  we  also  obtained  the  maximum  and  the 
average  theoretical  specific  capacity  of  Na-intercalated  M0S2  to  be 
146  mAh  g_1  and  86  mAh  g_1,  respectively.  Our  calculated  values  on 
the  specific  capacities  as  well  as  the  voltages  are  comparable  to  the 
recent  electrochemical  measurements,  which  reported  a  steady 
state  capacity  of  84  mAh  g_1  in  the  voltage  range  of  0.8— 1.2  V  [16], 
This  reasonably  good  agreement  between  the  measured  and 
calculated  electrochemical  characteristics  of  layered  MoS2  indicates 
an  ability  to  achieve  an  electrochemical  performance  close  to  its 
theoretical  maximum.  Furthermore,  our  results  on  electron  charge 
transfer  indicate  that  Na  atoms  donate  a  large  fraction  of  their 
valence  electron  to  2H  and  IT  polytypes  (see  Table  1).  This  large 
electron  donation  suggests  the  nature  of  the  bond  between  of  Na 
atoms  and  M0S2  to  be  of  the  ionic  character.  To  estimate  the  ease  of 


this  charge  transfer  during  cyclic  Na  insertion/deinsertion  in 
layered  M0S2,  we  have  used  chemical  hardness,  a  metric  that 
evaluates  cycling  performance  efficiency  [50].  Consistent  with 
earlier  calculations  [51  ],  a  low  value  (in  the  range  of  0.04—0.08  eV) 
for  the  chemical  hardness  of  the  Na  intercalation  in  layered  M0S2, 
point  out  to  the  easiness  of  charge  transfer  upon  Na  intercalation. 

One  of  the  crucial  characteristics  for  evaluating  suitability  of 
an  electrode  material  for  rechargeable  batteries  is  the  char¬ 
ge-discharge  rate,  which  depends  on  the  mobility  of  the  inter¬ 
calating  ion.  In  order  to  estimate  the  mobility  of  Na  in  M0S2,  we 
have  identified  its  diffusion  path  and  calculated  the  associated 
diffusion  barriers  using  the  climbing-image  nudged  elastic  band 
method  [52].  For  2H-  and  IT-M0S2,  we  considered  diffusion  of  an 
Na  ion  in  the  interlayer  van  der  Waals  spaces  for  concentrations 
x  =  0.125  and  x  =  0.5,  respectively.  Fig.  8c  shows  the  diffusion 
path  between  neighboring  low-energy  octahedral  sites  and  the 
associated  energy  profile.  In  both  cases,  the  path  passes  through 
the  high  energy  tetrahedral  binding  site,  with  diffusion  barrier 
being  0.68  eV  and  0.28  eV  for  2H  and  IT-M0S2,  respectively.  It  is 
evident  that  an  increased  Na— Na  interaction  at  higher  concen¬ 
trations  leads  to  a  faster  Na  diffusion  in  IT-M0S2.  The  Na  diffu¬ 
sion  barrier  in  IT-M0S2  is  smaller  than  Li  diffusion  barrier  in 
graphite  (0.4  eV)  [53]  and  is  comparable  to  Na  diffusion  barrier  in 
layered  cathode  materials  such  as  NaMPC^F  (M  =  Fe,  Feo.sMno.s, 
Mn)  and  layered  NaCoC>2  (0.3  eV)  [1,54],  Although  the  Na 
migration  barrier  in  2H-M0S2  is  larger,  2H-M0S2  still  offers  better 
mobility  in  comparison  with  a  number  of  electrode  materials 
such  as  post-spinel  phases  of  CaMn2C>4  and  CaT^CU  with  diffu¬ 
sion  barriers  in  excess  of  1.0  eV  [55],  Since  the  rates  of  diffusion 
in  layered  M0S2  is  limited  by  the  ability  of  van  der  Waals  spaces 
to  accommodate  Na  ions,  their  mobilities  can  be  further 
enhanced  by  using  composites  of  monolayer  M0S2  as  suggested 
by  recent  experiments  [20,21], 
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Fig.  7.  Electrode  potential  of  Na-intercalated  (a)  2H-MoS2,  (b)  1T-MoS2  against  Na/Na+. 


In  the  preceding  section,  we  have  focused  on  Na  intercalation  in 
layered  M0S2  in  the  concentration  range  of  0  <  x  <  1.  Recent  studies 
investigating  Li— M0S2  intercalates  have  reported  that  once  all 
interlayer  binding  sites  have  been  fully  occupied  with  Li,  further 
intercalation  at  higher  concentrations  (x  >  1)  leads  to  the 
destruction  of  layered  structure  of  M0S2  and  a  loss  of  cyclic 
reversibility  [33],  This  behavior  arises  due  to  inversion-type  re¬ 
actions  that  involve  breaking  of  Mo— S  bonds  and  the  formation  of 
lithium  sulphides  [33,15],  While  similar  inversion  reactions  have 
not  yet  been  reported  in  case  of  Na  intercalation  in  M0S2,  it  is 
nevertheless  worth  analyzing  their  reaction  energetics  considering 
close  magnitudes  of  the  electronegativities  of  Na  and  Li.  We 
therefore  calculated  the  reaction  energies  (A£r),  defined  as  the 
difference  between  the  energy  of  the  products  and  the  reactants,  of 
the  following  two  reactions:  Na  +  NaMoS2  ->  2NaS  +  Mo, 
A Er  =  -12.36  eV,  and  3Na  +  NaMoS2  ->  2Na2S  +  Mo, 
A Er  =  -15.46  eV.  Large  negative  energies  of  these  reactions, 
especially  when  compared  with  the  formation  energies  of  Na*MoS2 


Fig.  8.  Schematic  illustration  of  the  diffusion  path  for  Na  ion  in  (a)  2H-  and  (b)lT- 
MoS2,  from  a  low-energy  octahedral  (0)  site  to  its  neighboring  octahedral  site  via  a 
tetrahedral  (T)  site,  and  (c)  associated  energy  profiles  along  the  path  0  ->  T  ->  0. 


intercalation  compounds,  suggest  a  likelihood  of  Na-induced  bond 
breaking  in  the  host  M0S2  at  large  Na  concentrations  and  the 
subsequent  loss  of  charge— discharge  reversibility.  Severe  structural 
distortions  at  later  stages  of  Na  intercalation  in  M0S2  as  reported  in 
the  recent  experiments  could  possibly  be  attributed  to  such  Na- 
induced  inversion  reactions  discussed  here  [16],  In  light  of  the 
promising  electrochemical  properties  of  NaxMoS2  as  an  anode 
material  for  Na  concentrations  in  the  range  of  0  <  x  <  1,  these 
irreversible  reactions  highlight  the  need  for  limiting  Na  concen¬ 
trations  in  this  range,  and  motivate  the  development  of  inhibitor 
dopants  that  can  impede  the  onset  of  such  reactions. 


4.  Conclusions 

In  summary,  we  have  systematically  investigated  Na  intercala¬ 
tion  in  all  structural  polytypes  of  M0S2.  Na  ions  are  generally  found 
to  occupy  octahedral  or  distorted  octahedral  sites  within  the 
interlayer  van  der  Waals  spaces  upon  intercalation.  The  intercala¬ 
tion  of  Na  in  M0S2  causes  a  phase  transformation  at  x  =  0.39  from 
the  2H  to  IT  phase.  However,  we  find  that  Na  intercalation  in  M0S2 
is  thermodynamically  stable  across  the  entire  concentration  range 
of  0  <  x  <  1  against  the  phase  separation  of  Na.  Furthermore,  Nax_ 
MoS2  compounds  are  generally  metallic  with  a  low  average  elec¬ 
trode  potential  of  0.75—1.25  V  and  a  maximum  specific  capacity  of 
146  mAh  g_1.  Finally,  our  calculated  energy  barrier  of  0.68  and 
0.28  eV  for  Na  diffusion  in  2H-  and  IT-M0S2  suggests  a  moderately 
fast  charge  and  fast  discharge  rates,  respectively.  These  results 
suggest  that  IT  polytype  has  better  electrochemical  performance 
due  to  its  inherent  metallic  electronic  structure,  higher  capacity  for 
Na  ions  and  their  faster  mobilities. 

Since  the  capacity  of  layered  M0S2  is  limited  by  the  ability  of  van 
der  Waals  interlayer  space  to  accommodate  Na  ions,  exfoliating 
layered  structure  into  individual  monolayers  can  allow  for  an 
enhanced  adsorption  of  Na  ions  on  both  sides  and  faster  diffusion 
rates  as  hinted  in  some  recent  experiments  [20,21,56],  Similarly 
using  nanorstructures  of  M0S2  (such  as  nanosheets,  nanowires  and 
nanoribbons)  in  combination  with  composites  may  also  lead  to  an 
increase  in  Na  amount  adsorbed  per  M0S2  formula  unit  to  enhance 
the  specific  capacity.  Finally,  given  the  close  structural  and  chemical 
similarities  between  M0S2  and  other  transition  metal  dichalcoge- 
nides  such  as  WS2,  MoSe2  and  WSe2,  our  study  should  motivate  a 
thorough  investigation  of  Na  intercalation  in  several  transition 
metal  dichalcogenides.  Considering  that  graphite,  the  widely  used 
layered  material  as  electrodes  in  rechargeable  batteries,  is  not 
amenable  to  Na  intercalation,  our  study  suggests  the  family  of 
transition  metal  dichalcogenides  as  promising  candidates  for  an¬ 
odes  of  NIBs  for  stationary  energy  storage. 
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